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Electrons in a metal 



The Drude model was propose in 1900 by Paul Drude to 
explain the transport properties of electrons in a metal. 

This treats the electrons like an ideal gas. 
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It explains very well the DC and AC conductivity in metals, the 
Hall effect, and thermal conductivity. 

However, it greatly overestimates the heat capacities. 

http : / / en. wikipedia . org/wiki/Drude_model 
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Fermi-Dirac statistics 
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Why did the electron gas model fail? 



When we derived the Maxwell-Boltzmann distribution for the 
ideal gas, we assumed that it is extremely unlikely for two 
atoms to occupy the same energy level. 



The reason is that the energy levels are very close together, 
compared to the average energy of the atoms. 

Unfortunately, this is no longer true for electrons in a metal at 
room temperature. 
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Recall the equation for the energy of a state: 

h 2 2 2 2 

£n x n y n z = — ?( n x + n « + n z) 

Because the electron has a much smaller mass, the energy 
levels are much higher and further apart. 

So the electrons have to occupy all of the lowest levels first. 
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Fermi-Dirac statistics 



As a result, we cannot use the Maxwell-Boltzmann distribution 
for the ideal gas. 

Fortunately, we can use the same statistical methods that we 
have learnt so far. 

Lets start by looking at the picture of the bundle from another 
angle. 

. _. 



Q 

We have gi energy states in the bundle. Suppose the bundle 
contains n, electrons. 
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The exclusion principle 



Electrons are not allowed to occupy the same energy states. 
So they have to be stacked up from bottom to top. 
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When heated, most of the electrons are stuck - there is no 
space above to move up in energy I 

Only those few at the very top can. As a result, the heat 
capacities are much smaller than expected of a gas. 



Fermi-Dirac distribution 



According to the exclusion principle, exactly n, states must be 
filled and exactly g^ — states unfilled. 

The total number of ways the electrons can be arranged in this 
bundle is therefore: 

For the whole system - i.e. all the bundles - we get 



Applying the Lagrange multiplier method again, we would get 
rn 1 

exp(-Ai - X 2 si) + l 
where A x and A 2 are the Lagrange multipliers. 

We have applied the same constraints as before on the particle 

number N and the energy U. 
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Fermi-Dirac distribution 



Previously, we have used n(e) for the number density, g(e) for 
density of states, and de for the energy interval of the bundle. 

We also know that one of the Langrange multipliers would be 
related to temperature: 

A 2 = -1/fcflT 

The distribution of the electrons may then be written as 

n{e)de = ; ; r . 

exp(-Ai +e/k B T) + 1 

This is called the Fermi-Dirac distribution function. 
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Occupying the energy states 



For convenience, we define the following function: 

f(e) _ < £ ) 1 

<?(e) exp((-^ + e)/fc B T) + 1 

Since g(e)de is the number of states and n(e)de is the number 
of particles, /(e) would be the fraction of states that are 
occupied. 

So /(e) is called the occupation number. We need to get a feel 
as to what this looks like and how it changes with temperature. 

Lets start with the simplest case: T = OK. If we allow T to 
approach zero, we will find: 

/(e) = 1 for e < (i 
/(e) = 0 for e > fi 

This means that all states with energy below fi are fully 

occupied. All states with energy above ^ are empty. 
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Fermi-Dirac distribution 



We have obtained the Fermi-Dirac distribution function. 

niejde = . 

exp(-A 1 +e/fc i3 T) + l 

The total number of electrons is fixed and is given by 

integrating the normalisation condition: 

N _ r°° g(e)tfe 

Jo exp(-A x + e/k B T) + 1 
In principle, we could solve for the Lagrange multiplier Ax- It is 
common practice to express it in terms of fj, as follows: 

exp(-Ax) = exp(~^/k B T) 
H is called the "chemical potential." 



Occupying the energy states 



The graph for /(e) looks like this 
f(s) 



The shape shows that the energy levels are occupied below a 
certain energy, and unoccupied above that. 



This feature is characteristic of the Fermi-Dirac distribution 
that we are studying. 
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The Fermi energy 



At T = OK, The highest energy in the stack of electrons is fx. 
This energy is also called the Fermi energy, Ep. 



We can find the Fermi energy E F by integrating the 
normalisation condition: 



N = 2 x 



roo 

/ n(e)de 
JO 



The factor of 2 must be added because each energy level can 
be occupied by 2 electrons - spin up and spin down. 

We shall solve this for the Fermi energy Ep. 



The Fermi energy 

We have previously derived the density of states: 

4-mirV , 



9(e) 



-(2me) 1 / 2 



In the topic on ideal gas, this is obtained by counting the 
number of energy states of a particle in a 3-D box. In this 
topic on electrons, we have used the same particle in a box. 
The difference from ideal gas only arises later on, when we 
make different assumptions about the energy levels. 

So the same formula for the density of states can be used for 
both the ideal gas and the electrons. We can therefore 
substitute the formula into the normalisation integral 

rE F 

N = 2 x / g(e)de 
and solve for the Fermi energy. The result is 

h 2 {3n 2 N\ 2/3 



Ej, = 



2m 
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The Fermi energy 



In terms of the occupation number, 

n(e) = g(e)f(e) 
So the normalisation condition can be written as: 

N = 2 x J n(e)de = 2 x J g(e)f(e)de 

We know that at OK, /(e) = 0 for e > fi. So So the integration 
would stop at e = fi: 

N = 2 x J^ F g(e)f(e)de 

since (j, = E F at OK. 

We also know that at OK, /(e) = 1 for e < fi. So 

rE F 

N = 2 x / g(e)de 
Jo 
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Electronic heat capacity 
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Heat capacity 



We want to use the ideas and formulae that we have developed 
to calculate the heat capacity of electrons. 

To see how to do this, recall that the Drude model would 
predict a heat capacity for that is the same as that of an ideal 
gas: 

C = -Nk B . 
2 B 

It is known from experiments that the actual heat capacity of 
the electrons is much smaller. This refers to measurements 
that are done at room temperature. 
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Heat capacity 



We start by assuming that the thermal energy is indeed much 
smaller than the Fermi energy. We shall derive an expression 
for this thermal energy, and then calculate it at room 
temperature to see if the assumption can be justified. 

At a temperature above 0 K, the occupation number 



exp((- M + e)/fc B r) + 1 
would no longer have a sharp step at the Fermi energy. 




If k B T is much smaller than fi, the graph would remain close to 
the step, as if the step has become smoother. 
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Heat capacity 



We can understand this if we think that the electrons that are 
stacked up to the Fermi energy do not have enough energy to 
jump out of the stack. 



C 
LU 



This would only be true if the thermal energy is much smaller 
than the Fermi energy. 

In order to find out if this is true, we need to estimate this 
thermal energy. 
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Heat capacity 

The "smoothened" slope of the graph tells us that electrons 
just below the Fermi energy (/x = E F ) is excited just above it. 



So we can estimate gain in thermal energy of the excited 
electrons by the width 5 of this slope. 



In order to do this, we take a closer look at the occupation 
number 

f(e) 



exp[(e- M )/fc B T] + 1 
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Heat capacity 



We are assuming that k B T is much smaller than p. 

For energy e higher than /i by a few times of k B T, The 
exponential function exp[(e — n)/k B T] would quickly become 
large. 



"he occupation number is then approximately 



exp[(e - n)/k B T] + 1 exp[(e- M )/fc B T] 
which is just the exponential function with negative argument 

/0) = exp[-( e -/x)/fc B T]. 

m exp[-(E-n)/k B T] 
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Heat capacity 



There is a corresponding behaviour just below the Fermi energy. 

For energy e lower than fj, by a few times of k B T, The 
exponential function exp[(e — ^)/k B T] would quickly become 
small. 

The occupation number is then approximately 

exp[(e - n)/k B T\ + 1 
where we have used the binomial expansion and kept only the 
first order term. 

f(f) 1-exp[(E-m/k B T| 

I 
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Heat capacity 



This means that the part of the graph to the right of fi falls off 
exponentially. 



f(s) 




It falls by a fraction of 1/e over an energy range of k B T. 
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Heat capacity 



This means that the part of the graph to the left of fi tends to 
the line /(e) = 1 exponentially. 




It reaches within /(e) = 1 by a fraction of 1/e, over an energy 
range of k B T. 

This means that it is the electrons within this energy range 
that is excited. So the thermal energy of the excited electrons 

is of the order of fc R T. 
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Heat capacity 



At this point, we should justify our assumption that k B T at 
room temerature is much smaller than the Fermi energy 
(i = E F , which is defined at 0 K. 

We shall take sodium metal as an example, and calculate k B T 
and E F for this metal. 

In sodium, each atom has one valence electron. This electron 
is mobile and forms the electron gas that we are talking about. 
In order to calculate the Fermi energy, we need the number 
density N/V. We can calculate this from the following data: 

density = 0.97 g cm -3 

relative atomic mass = 23.0 

So the volume for one mole of atoms is 

23 -f- 0.97 = 23.71 cm 3 , 
24 



Heat capacity 



In contrast, at room temperature 298 K, we can calculate that 

1 

k B T w — eV. 
u 40 

This is about 120 times smaller than the Fermi energy. 



We can repeat this for other typical metals, and we would get 
similar answers. 

This justifies our assumption that at room temperature k B T is 
much smaller than the Fermi energy. 

We are now a step closer to estimating the electronic heat 
capacity. Next, we need to understand the behaviour of the 
excited electrons. 
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Heat capacity 



and the number density is 

N 

— = N A + 23.71 

V A 

where N A is Avogadro's number. This gives an answer of 
2.54 x 10 28 m" 3 . 



Using the Fermi energy formula 

Ep 



h 2 /3^ 2 A^ 2/3 



2m \ V 

where m is the mass of the electron, we can find that the 
Fermi energy is 3.16 eV. 
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Heat capacity 



When temperature increases above 0 K, the step in the 
Fermi-Dirac distribution becomes smoother as electrons just 
below the Fermi energy are excited above it. 




Notice that the "tail end" of the distribution - to the right - 
looks exponential. This is because there are relatively few 
electrons above the Fermi energy. So these can behave like the 
ideal gas and approximately obey the Maxwell-Boltzmann 
distribution. 
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Heat capacity 



It is clear from the graph that for a small increase in 
temperature, only electrons close to the Fermi energy are 
excited. Most of the electrons are below the Fermi energy and 
are not excited at all. Since they cannot be excited, these 
electrons would not contribute to the heat capacity . 




So it is mainly the electrons close to the Fermi energy that 
would contribute to the heat capacity. So we can use these to 

estimate the heat capacity and ignore the rest. 
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Heat capacity 



The electrons that do get excited are in the small energy range 
of order k B T from the Fermi energy E F . 




So we can estimate N\ with the number of electrons that are 
within the energy interval of 

de = k B T 

from the Fermi energy. 
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Heat capacity 



As we have seen, these electrons close to the Fermi energy 
behave like the ideal gas. We know that the energy of an ideal 
gas is 

where N± is the number of particles in the ideal gas. 

In the case of the electrons, iVi should refer to the number of 
electrons above the Fermi energy, and not the total number of 
electrons. We can estimate this number as follows. 
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Heat capacity 




We know that the number of particles in a given energy 
interval is 

n(e)de = 2g(e)f(e)de 

where the factor of 2 again comes from the spin states of the 
electrons. 

At OK, the energy states below E F are fully occupied, i.e. 
/(e) = 1. So the number would given by 

iV x « 2g(E F )k B T 
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Heat capacity 



From the definitions of the density of state and the Fermi 
energy, it can be shown (see exercises) that 

g(E F ) 



4E F 

where N is the total number of electrons. This gives 

\AEpJ B 2 Ep 
Substituting into the energy for ideal gas 

Ui = |jVifc B T 

we get the energy for the excited electrons 

3 (3 k R T\ 9 k. n r 2 

Ui = -[ -N-^— k B T = -Nk 



2 \2 E F J 4 E F 

Differentiating with respect to T, we get the electronic heat 
capacity 

C = 9 -Nk B k -^ 

2 _ E F 
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Heat capacity 

1. Only electrons within k B T of the Fermi level can jump up. 
Those further down cannot jump because the states above are 
still mostly occupied. 

2. The number of electrons that can get excited is ^pjyN. 
This is a fraction ^ fc jJ of the total. 

3. Each would gain an energy of about 3k B T/2. So the 
increase in total energy is obtained by multiplying the number 
of excited electrons by this energy 

AU « x 3k B T/2 = 9 -Nk B k -^. 

2E F B 1 4 U E F 

4. Differentiating with respect to T, we get the electronic heat 
capacity 

9 k R T 

C = -Nkg-*— 

2 u E F 
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Heat capacity 



We summarise what we have learnt about the heat capacity: 

Some states below Fermi level becomes empty, and some 
states above becomes occupied. 



For a temperature T, an electron that gets excited would be 
able to gain on average an energy of about 3k B T/2. 



Heat capacity of a metal 



We have obtained the electronic heat capacity 

C = 9 Nk B k -^ 
2 U E F 

More detailed calculations show that the factor of 9/2 should 
really be ?r 2 /2: 

C = *-Nk B k J^ 
2 u E F 

This equation is often written in the form 

7V 2 T 

C = —Nk B — 
2 U T F 

where T F = E F /k B is called the Fermi temperature. 
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Heat capacity of a metal 



Notice that the heat capacity is directly proportional to T. 

7V 2 T 

C = —Nk B — 
2 Tp 

This is often written in the form 

C = jT. 

There is another contribution to the heat capacity. This comes 
from the vibrations of the atoms, and it is proportional to T 3 . 

We could imagine writing the total heat capacity in the form: 

c v = -1T + AT 3 . 

We can measure this heat capacity to check if the formula is 
correct. Suppose that we have obtained a table of values for T 
and cy. To check if the formula is correct, we can rewrite it in 
this form: 

cy/T = 1 + AT 2 

If we plot cy/T against T 2 , we should get a straight line. 
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Effective mass, and a few other things 
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Heat capacity of a metal 



In 1955, William Corak and his fellow co-workers measured the 
heat capacities of copper, gold and silver from IK to 5K in 
their laboratory in Pittsburgh. 

They plotted cy/T against T 2 , and got the straight lines. The 
picture here is a sketch of their results. This shows that the 
predictions of the Fermi-Dirac statistics are correct. 



Heat capacity of silver 




5 10 



T 2 (K 2 ) 
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Heat capacity of a metal 



From the formula, 

cv/T = 7 + AT 2 

we know that we can find 7 from the y-intercept of the graph. 
The following table shows the values for a few metals. 



FREE ELECTRON y MEASURED >' RATIO" 

element (in 10" 4 cal-mole" J -K" 2 ) (m*fm) 

Li 1.8 4.2 2.3 

Na 2.6 3.5 1.3 

K 4.0 4.7 1.2 



The second column contains the values of 7 from the formula: 

C = ^iVfc B fc ^ = 7T 
2 u E F 

The third column contains the values actually measured. They 

are obviously different. 
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Heat capacity of a metal 



Does this mean that our theory is wrong? 

Yes and no. The theory tells us that the electronic heat 
capacity is proportional to T. The measurement shows that 
there is indeed such a contribution. 

So we are not completely wrong. Perhaps the theory needs 
refining. We can be optimistic and go back and try and 
understand what we have missed. 

Recall that we have started with a particle in a 3-D box and 
calculated the energy levels. Then we just fill these up with 
electrons and calculated the heat capacity. All we have is a gas 
of electrons in empty space. 
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Heat capacity of a metal 

If we imagine that the electrons in a metal has a different 
effective mass m* than its natural mass, we can "explain" the 
difference in 7. 

According to the formula for heat capacity, calculated 7 is 
given by 

This is inversely proportional to the Fermi energy 

_ h 2 f3n 2 N\ 2/3 

which is in turn inversely proportional to the mass m of the 
electron. 

So 7 is directly proportional to the mass. 
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Heat capacity of a metal 



But what about the atoms? The "empty space" is really filled 
with atoms. The electrons must surely interact with the atoms. 
This then is the reason for the difference between theory and 
measured 7. 

However, the measured C is proportional to T. This agrees 
with theory, and it should mean something. One possibility is 
that, for some reason, the electron interacts only weakly with 
the atoms. This idea has been shown to be correct by other 
types of measurements. 

According to this idea, the behaviour of the electrons in the 
presence of the atoms is essentially the same. The difference is 
that the interaction with the atoms make the electrons behave 
as if they have a different mass. 
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Heat capacity of a metal 



This means that we can get the effective mass if we divide the 
measured 7 (in the third column) by the calculated 7 (in the 
second column). 



This is shown in the last column of the table. 



FREE ELECTRON y MEASURED f RATIO" 

element (in 10~ 4 cal-mole~ ^K" 2 ) (m*/w) 

Li 1.8 4,2 2.3 

Na 2.6 3.5 1.3 

K 4.0 4.7 1.2 
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Behaviour at high temperatures 



At very high temperatures, all the electrons could get excited. 
Then they would start to behave like an ideal gas. There are 
many more energy state they can reach and they are much less 
likely to be forced towards the same energy states. 

This would happen only if the electrons are excited far above 
the Fermi energy. Therefore we can use Fermi temperature as 
a reference. A temperature would be high if it is high compared 
to the Fermi temperature. 

At the high temperature, the heat capacity would therefore 
change to that of the ideal gas. So instead of 

TV 2 T 

C = —Nk B — 
2 a T F 

which is very small, it would become 

C = -Nk B 
2 a 

which is much larger. 
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A note on fermions 



Note that electrons belong to a larger family of particles called 
fermions. 

A fermion is a particle with a half integer spin. Other examples 
are 

the proton, 
the neutron, 

the helium-3 ( 3 He) atom and 
the oxygen-13 ( 13 0) atom. 

All fermions are known to show the kind of properties we have 
seen for electrons. They obey the Fermi-Dirac statistics. 
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A note on the chemical potential 



We have seen many times that the chemical potential ^ is 
equated to the Fermi energy E F . It is important to note that 
this is true only at 0 K. 

Recall the chemical potential fi is determined by the number of 
particles: 

Jo exp[(e - n)/k B T] + 1 
At T = 0 K, the integral is just 

N = j^g{e)de 

which is why /x is the Fermi energy then. 

At higher temperature, fj, would change. But for temperatures 
well below the Fermi temperature, as with electrons in metals 
at room temperature, we may assume that fi is still quite close 

to the Fermi energy. 
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What we have learnt so far 

1. Electrons follow the Fermi-Dirac distribution: 

, s , g(e)de 

n(e)de = — — . 

exp[(e - n)/k B T] + 1 

2. At 0 K, the electrons fill up the energy states from the 
lowest level. The highest energy of the electrons is called the 
Fermi energy: 

h 2 {3n 2 N\ 2/3 

E F = 

2m\ V 

3. When thermal energy k B T is small compared to the Fermi 
energy, only electrons near the Fermi energy would get excited. 
So only these would contribute to the electronic heat capacity: 

n 2 T 
C = —Nk B — 
2 *T F 

4. The Fermi temperature is defined as 

k B 
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Worked Examples 
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Example 

We can make use of the following relations for the density of 
states. 

Relation with the total number of states: 

as 

Formula for the total number of states: 
so 

GO) cx e 3 ' 2 . 
Relation with the Fermi energy and N: 

N = 2 j g(e)de. 
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Example 



Example 1 

Show that the density of states at the Fermi energy E F is given 
by 

, « 3iV 

g(E F ) 



where N is the number of electrons. 



From 



and 



we have 



From 



we have 
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Example 



, x G(s) 
els 



N = 2 g(e)de, 

Jo 

N = 2G(E F ). 

G(e) oc e 3 / 2 , 
9(e) _ 3 



G(e) 2e 
At the Fermi energy E F , we get ... 
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Example 



9{Ef) ~ ^EfT- 
Substituting the above result for N: 

N = 2G(E F ), 

we get 

^ 3iV 

g(E F ) = . 

AE F 
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Example 

We need to maximise the Lagrange function In Q + X\N + A 2 £/. 

Using Stirling's formula 

In AT! ks N\nN — N 

we have 

In Q w In ffi-Sj-(rij Inrii-nJ - [(ffj-nj) In^-n,) - (gj-rij)] 

i 

Differentiating with respect to rij, 

9 In Q — tii 
— = - In m + InCffi - nO = In - 

on-i n i 

The total number is given by 

i 

and the total energy by 

u = Yl n i e i- 

i 

Differentiating these ... 
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Example 



Example 2 

The number of ways of arranging fermions is given by 



Use this formula, along with the conservation of particle 
number and total energy, to derive the Fermi-Dirac distribution 
function: 

rh 1 

gi exp(-Ai - A 2 e,) + 1 
where A x and A 2 are the Lagrange multipliers. 
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Example 



We get 

dlniV 



= 1 

dm 



and 

d\nU 



drii 

Then differentiating the Lagrange function, we get 

Rearranging, we get the Fermi-Dirac distribution function: 

n,- 1 



gi exp(-Ai - A 2 £i) + 1 
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Example 



Example 3 

Show that the mean energy of electrons at absolute zero, e is 
3E F /5, where E F is the Fermi energy. 
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Example 



Integrating, 



To relate this to N, 



rEp 

N = 2 g(e)de. 



we substitute the density of states 



and integrate: 



N = 2 =—(2meY' 2 de 

Jo h 3 K ' 



N = 2^(2^)3/^4/2 



The average energy is the energy per particle 

- U 3^ 

e = — = —Ep 
N 5 



58 



Exercises 



The total energy is 

TOO 

/ n(e)ede. 
JO 



/•oo 

U ■ 



For the Fermi-Dirac distribution, 



U = 2 i g(e)f(e)ede. 
At 0 K, this would be 



JO 



r E F 

U = 2 g(e)ede. 



Recall the density of states 



g(e) = — fer (2m £ ) 1 / 2 . 



Substituting into the integral, 



U = 2 ^—(2me) 3 > 2 de. 

JO h-3 
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Example 



Example 4 

Metallic hydrogen can be produced by compressing liquid 
hydrogen to extremely high pressure. This transition has been 
observed when the density reaches 0.32 x 10 6 mole m~ 3 , at a 
temperature of 3000 K. It is thought that in the metallic phase 
the liquid comprises Hj" ions and electrons. Were the electrons 
in this liquid metal obeying classical or Fermi-Dirac statistics. 
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Example 



We need to compare the 3000 K temperature to the Fermi 
temperature Tp. 

If the temperature is much higher than Tp, the electrons would 
obey classical statistics and behave as an ideal gas. If the 
temperature is much smaller than Tp, the electrons would obey 
Fermi-Dirac statistics. 



"he Fermi temperature is related to the Fermi energy by 

E F 
k B 



and the Fermi energy is given by 

Ep 



r? /3tt 2 aa 2/3 



2m \ V 

where m is the mass of the electron. 
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Example 



Example 5 

One possible means of obtaining fusion energy is to implode 
spherical capsules containing heavy hydrogen by irradiating 
them with high power lasers. For fusion to occur the implosion 
core needs to reach electron temperatures of the order of ~ 1 
eV, and the electron density needs to be of the order of 
10 33 m -3 . Do we need to use Fermi-Dirac statistics to describe 
the electrons, or are Maxwell-Boltzmann statistics sufficient? 
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Example 



The hydrogen density 0.32 x 10 6 mole m -3 is also the electron 
density, since each molecule H 2 converted to Hj" gives one 
electron. The number density is obtained by multiplying 
Avogadro's constant: 

— = (0.32 x 10 5 ) x (6.022 x 10 23 ). 

Substituting this and the other constants into the above 
formulae, we get the Fermi temperature 

T F = 141,0007f. 

The temperature of 3000 K is therefore much lower than the 
Fermi temperature. So the electrons in metallic hydrogen 
would obey Fermi-Dirac statistics, in the same way as normal 
metals at room temperature. 
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Example 



We need to compare the 1 eV electron temperature to the 
Fermi energy E F . The Fermi energy is given by 

t 1m\ V ) 

The number density N/V is 10 33 m -3 . Substituting into the 
above formula, we get 

E F = 3650 eV. 

The electron temperature of 1 eV is much smaller than this. 

Therefore the electrons would be described by Fermi-Dirac 
statistics. They would behave like electrons in metal at room 
temperature, and not as an ideal gas. 
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Example 



Note 

Notice that the terms energy and temperature have been used 
interchangeably. 

A temperature of E = 1 eV, means that the temperature 
would be 

T = E/k B 

in Kelvin. 
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